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Small distortions in the observed shapes of distant galaxies, a cosmic shear due to gravitational 
lensing, can be nsed to simultaneously determine the distance-redshift relation, r{z), and the density 
contrast growth factor, g[z). Both of these functions are sensitive probes of the acceleration. Their 
simultaneous determination allows for a consistency test and provides sensitivity to physics beyond 
the standard dark energy paradigm. 


Introduction. The observed acceleration of the cosmo¬ 
logical expansion is driving a revolution in fundamental 
physics. This revolution could transform our understand¬ 
ing of particles and fields (through the discovery of a new 
ingredient, the “dark energy”) or revise our deepest un¬ 
derstanding of space and time (by forcing fundamental 
changes to our theory of gravity). In this Letter we dis¬ 
cuss how wide and deep tomographic cosmic shear sur¬ 
veys, through their sensitivity to both geometry and the 
growth of density perturbations, can be used to distin¬ 
guish between these two possibilities. We also emphasize 
the broader utility of having these two probes of, or win¬ 
dows on, acceleration and gravitation. 

Despite the variety of phenomena that can be ex¬ 
plained with the cold dark matter cosmolo^ augmented 
with a dark energy component HiiilEli , we still 
only know of dark energy through its gravitational in¬ 
fluence. And unlike dark matter, we have little hope of 
directly detecting the dark energy via earth-bound labo¬ 
ratory experiments. 

Given that we only know of dark energy through its 
gravitational effects, we must bear in mind the possibility 
that what we explain with dark energy, may actually be 
due to corrections to Einstein gravity. Note, as a histor¬ 
ical precedent, that the anomalous perihelion precession 
of Mercury detected in the 19th century was first ex¬ 
plained with unseen matter 0 before Einstein provided 
the correct explanation. 

Assuming Einstein gravity, the growth of cold dark 
matter density contrasts in the linear regime ((5(x, t) = 
Sp{x,t)/p{t) « 1 ) can be written as S{x,t) = 

g{t)6{x,ti) where U is some early time and g{t) is called 
the growth factor, usually written as a function of red- 
shift, z, instead. The growth of density contrasts re¬ 
sults from a competition between the gravitational force 
pulling matter toward overdensities and the expansion of 
the Universe driving everything apart. Thus g{z) is sen¬ 
sitive to both the gravitational force law and the history 
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of the expansion rate. With the history of the expansion 
rate determined by r[z), g{z) can then be used to test 
the gravitational force law on Mpc and larger scales and 
thereby distinguish Einstein gravity from alternatives. 

More generally, inconsistency of the standard dark en¬ 
ergy paradigm with the combination of r(z) and g{z) 
could arise for a variety of reasons. For example, the 
growth factor could also be altered by non-gravitational 
interactions of the dark matter. The cosmic-shear in¬ 
ferred r(z) and g{z) may be internally consistent, but 
inconsistent with tTz) as inferred from supernovae due 
to axion-dimming |^. It is thus imperative to probe ge¬ 
ometry and growth in as many ways as possible. 

Cosmic Shear Basics. Weak gravitational lensing maps 
source galaxies to new positions on the sky, systemati¬ 
cally distorting their images. The resulting shear 7 of 
their images is related to the projected foreground mass 
contrast inside an angular radius 9: 

7 t = k{< 0) — k{9) where 7 * is the tangential component 
of the shear, k is the mass surface density divided by 
(c^/47rG)(rs/rins) and the are the angular diameter 
distances of the source, lens, and lens-source 0 . 

By separating the galaxies into n redshift bins, la¬ 
beled by i, we can create n shear maps, 7 ^. The most 
interesting statistical property of these maps, and the 
sole one we will consider here, is the two-point func¬ 
tion, (jijj). This two-point function is most easily ex¬ 
pressed in the spherical harmonic space in which we have 
(7^7j™) = Cl^SivS^m'. 

These n{n + l)/2 unique shear power spectra can be 
written as a projection of the matter power spectrum, 
A'^{k,z) along the line of sight: 

Cf = 'k‘^1/2 J drrW{fi,r)W{fj,r)A'^{k,z{r)) ( 1 ) 

where IU(f,-,r) = for r < r, and zero otherwise 

|inl lll| . Here t; is the distance to galaxies in redshift 
bin i. The power spectrum is evaluated at the redshift 
corresponding to distance r on our past light cone and at 
wavenumber given by fc = l/r. For simplicity we approx¬ 
imate the distribution of galaxies in each redshift bin as 
a spike of zero width at the center of the bin. Below we 
assume n = 8 redshift bins centered on Zi = 0.2 -|- 0.4i 
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where i runs from 0 to 7. In the top panel of Fig. 1 we 
show the auto power spectrum for sources at z = 1. 
Reconstructing r{z) and g{z): Qualitative. With cosmic 
shear maps from multiple source redshift bins, one can 
simultaneously determine both g{z) and the distance- 
redshift relation, r{z). Counting the degrees of freedom 
one can see that the multiple source redshift bins are 
essential: it would be impossible to reconstruct two free 
functions from just a single shear power spectrum. Using 
multiple source redshift bins provides us with the neces¬ 
sary further constraints. 

To gain further insight into the reconstruction, con¬ 
sider that the shear power spectrum of a given source 
redshift bin is a sum of shear power from lenses over 
a range of redshifts, as illustrated in the top panel of 
Fig. 1. Increasing g(z) simply increases the amplitude of 
the shear power contribution from structures at redshift 
z by g‘^{z). Increasing r{z) also changes the amplitude 
of the contribution from lenses at redshift z. In addition 
it causes a shift of the power toward higher I since the 
three-dimensional structures, if further away, will project 
into smaller angular scales. For a single source redshift, 
the changes in the shear power spectrum due to a change 
in r(z) at just one redshift, could also come from the 
appropriately chosen changes to g{z) over a range in red- 
shifts. Thus it is impossible to simultaneously recon¬ 
struct g{z) and r(z) from the shear power spectrum of 
a single source redshift bin. Including multiple source 
redshift bins breaks this degeneracy. 

If the matter power spectrum were a power law (and 
therefore all the curves in Fig. 1 were power laws), then 
we would still have a degeneracy between growth and 
distance even with multiple source redshift bins. The si¬ 
multaneous determination is enabled by a bend in the 
matter power spectrum, the exact location of which de¬ 
pends on the matter density today, pm ■ This feature, 
calibrated by CMB determination of pm, acts as a ‘stan¬ 
dard ruler’ [l^. 

Reconstructing r(z) and g{z): Quantitative. Our proce¬ 
dure is a modification of that in [llll where simultaneous 
reconstruction of distances and growth factors from cos¬ 
mic shear data was first considered. 

We parameterize r(z) by its values specified at discrete 
redshift values Zi = 0.4i for i = 1 to 8. In addition we 
set r(z*) = Vs/0s where z* ~ 1089 is the redshift of 
last-scattering. The sound horizon at last-scattering, Ts, 
depends on pm and pb. These and the angular size of the 
sound horizon, Og, are constrained by CMB observations. 
The values of r at all other points of z are found by linear 
interpolation. The A^(fc, z(r)) factor in Eq. n]is evaluated 
by inverting r(z) to get z(r). 

We assume that the primordial curvature power spec¬ 
trum, with amplitude at wavenumber k specified at hori¬ 
zon crossing (when k/a = H), is a quasi-power law 
with logarithmically varying spectral index, nsik) = 
ns + cts\n{k/kf) where kf = 0.05Mpc“^. This is the 
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FIG. 1: Dependence of the z = 1 shear auto power spectrum 
on g{z) and r{z). There are n — 1 more auto power spectra 
and n(n — l)/2 cross power spectra not shown here. The 
solid line is the shear power spectrum for sources at z = 1. 
The dashed lines show the contributions to this shear power 
spectrum from lens slices of width Az = 0.2 centered at z = 
0.1, 0.3, 0.5, 0.7 and 0.9. Their sum gives the solid line. The 
lower panel shows the z = 0.5 contribution again (dashed 
line), how it would look with an increase in g{z = 0.5) (dotted 
line) and how it would look with an increase in r(z = 0.5) 
(dot-dashed line). 


form of the expected power spectrum from inflation. The 
power spectrum at fixed time (or redshift) is related 
to this primordial power spectrum by a scale-dependent 
transfer function, T(fc) and a growth factor, g{z) so that 

P7r2 

AL(fc,z) = -^A,{k/kfr^^^)-^g\z)T\k) (2) 

In general the time and scale-dependence do not factor as 
written here, but we are interested in sufficiently small 
scales where the dark energy perturbations can be ig¬ 
nored and in this case all modes grow at the same rate. 
The ‘lin’ subscript on here stands for ‘linear theory’. 
We take non-linear evolution into account using the pre¬ 
scription of Peacock and Dodds . 

We parameterize g"^{z)/a? {z)^^ by its value at the 
eight discrete redshifts used for parameterizing r(z), plus 
its value at z = 0. We assume g(z*) = a(z*), as is the case 
in the matter-dominated era, and calculate 5 ^(z)/a^(z) 
at non-grid values of z by linear interpolation. 

The transfer function above depends on the matter 
content. In addition to Aq^ ns and as, the shear power 
spectra are also therefore affected by pm, Pb and the en- 
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ergy density of the cosmological neutrino background. To 
control these contributions we assume we have a measure¬ 
ment of the CMB temperature and polarization power 
spectra as expected from Planck. Since these CMB power 
spectra are also affected by the redshift of reionization, 
Zri, and the primordial fraction of baryonic mass in "^He, 
Yp, we include these parameters as well. Our parameter 
set is thus 0s, Aq, ns, as, Pm, Pbi ^rij Yp), eight 
r(z) parameters and nine g^{z)/a?'{z) parameters. 

To forecast errors, we Taylor expand to first order the 
dependence of the shear power spectra on these param¬ 
eters about our fiducial model. The expected covariance 
matrix for the errors in the estimated parameters can 
then be calculated via, e.g., Eq. 21 of 
Acceleration Without Dark Energy. As an example of 
acceleration without dark energy we turn to the ‘self- 
inflating’ branch of the DGP model 0 In this model 
our 3-|-1-dimensional world (or ‘brane’) is embedded in a 
4-|-1-dimensional space. The Friedmann equation on the 
brane becomes 

- H/r, = ^Pm (3) 

and thus H tends to a constant (1/rc) as the Universe 
expands, just as it would with the usual Friedmann equa¬ 
tion in the presence of a cosmological constant. 

The extra dimension also leads to a modification of 
the Poisson term on scales between a smaller scale that 
is perhaps about 1 Mpc and a large scale, rc- Song 0 
shows this to be 
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where S = Spm/pm [13 ■ Growth is suppressed in DGP 
gravity relative to dark matter only Einstein gravity by 
the factor in square brackets. However, dark energy also 
suppresses growth by Pm/ptot, which, for Einstein grav¬ 
ity, replaces the factor in square brackets. In the fol¬ 
lowing we set Tc = 1.27/i?oll3- For non-linear growth 
we apply the Peacock and Dodds prescription, as we do 
for Einstein gravity, although consequences of the DGP 
model for non-linear growth are at the moment unclear. 
Quantitative Forecasts for a Fiducial Survey. We take 
the fiducial survey to be the “27r” deep wide survey of 
20,000 square degrees in six wavelength bands from 0.4- 
1.1 /im to be undertaken by the Large Synoptic Survey 
Telescope (FSST). Several hundred sky-noise limited ex¬ 
posures in each optical band will be obtained for each 10 
square degree sky patch over a period of ten years. The 
shapes of galaxies out to a redshift of 3 (integrated galaxy 
density of 50 per square arcminute) will be measured at a 
precision far exceeding the 0.15 intrinsic random shear of 
an individual source galaxy. These galaxy redshifts will 
be estimated from fits of the 6-band fluxes to spectral 
templates for galaxies vs type and redshift, a technique 
called photometric redshift estimation. 


This survey will yield the shear and redshift of 3 billion 
source galaxies over a redshift range of 0.2 - 3. Based on 
experiments with the active optics 8m Subaru telescope, 
systematic shear error will be kept below 0.0001 on all 
angular scales considered here. 

Here we model the noise in the resulting shear maps as 
in 0 Our analysis does not include effects of redshift 
errors. For systematic errors in distance determinations 
to be less than 1%, it is sufficient to require Az/z < 0.01 
(since r varies slower than linearly with z) where Az is 
the error in the mean redshift of a given redshift bin. 
Simulations based on current surveys indicate that this 
level of control is achievable at z > 0.1 0 

Galculating sufficiently accurately at small scales 
is difficult [ 3 . We conservatively discarded data 
at I > 1000. We have not modeled fluctuations in the 
dark energy, which can be important on larg e scales, and 
therefore discard data with ^ < 40 as in |l5|. 

Results. Our results are presented in Fig. 2. To simu¬ 
late reconstructed r{z) and g{z) we add a realization of 
the errors, drawn from a zero mean multivariate Gaus¬ 
sian with our forecasted covariance matrix, to the fiducial 
values for r(z) and g{z). The error bars are the square 
root of the diagonal elements of this covariance matrix. 

Distances are reconstructed with ~ 2% errors, even 
out to z = 3.2. The distance errors are highly correlated; 
certain linear combinations will have even smaller errors. 
The growth factors have 3% to 4% errors at z < 1.2 and 
then grow steadily with z. The tightness of the constraint 
at z = 3.2 is an artifact of our parameterization. It is due 
to the g^(3.2)/a^(3.2) parameter influencing g{z) all the 
way out to z* because of our interpolation scheme. 

One can distinguish dark energy from the DGP model 
even if the dark energy density evolution is adjusted to 
match r(z) for the DGP model. As mentioned above, 
these two scenarios will make different predictions for 
g{z). In the right panels the g{z) curves for the fiducial 
DGP model and for the dark energy model with identical 
H{z) are shown. Their difference in values is 221, 
corresponding to almost 15cr. 

As a test of our calculations, we used our covariance 
matrix for r(z) and g{z) to calculate constraints on w, 
assuming a dark energy model with constant w = P/p. 
The results agree with a more direct calculation of the 
expected error in w that bypasses the g and r param¬ 
eterization. The constraint on w is almost entirely due 
to r constraints rather than g constraints. The roles of 
geometry and growth are also discussed in [i^, 0| . 

We note that these measurements of distances into the 
matter-dominated era, combined with Planck’s GMB ob¬ 
servations can be used to achieve (T((Dtot — 1)^^) 10”^ 

|22|. greatly improving the precision with which this ro¬ 
bust prediction of inflation can be tested. Allowing for 
non-zero curvature will mean just one more parameter to 
fit in our analysis and so will not qualitatively degrade 
our eight distance determinations. 
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FIG. 2: Reconstructed distances (left panels), and growth 
factors (right panels). The lower left panel shows the frac¬ 
tional residnal distances, \r{z) — where r(z) 

are the reconstrncted distances and rfld(a) are the distances 
in the hducial DGP model. The lower right panel shows the 
residual growth factor, g(z) — giid(z). The curves in the right 
panels are gndiz) (solid) and g(z) for the Einstein gravity 
model (dashed) with the same H{z) and pm as the DGP 
model. Although these two models have the same r{z) they 
are distinguishable by their significantly different growth fac¬ 
tors. 


Discussion and Conclusions. The statistical errors in our 
fiducial survey are small enough to allow very precise 
reconstructions of distance and growth as a function of 
redshift. We have argued that redshift errors will not 
qualitatively affect our results. 

To investigate the impact of shear calibration errors, 
we parameterized the observed as fifj times the true 
C\^ and extended our parameter set to include one gain 
parameter, fi, for each source redshift bin, i. With a 
1 % prior determination of all the calibration parameters, 
we find that the distance errors increase by less than 
25%, growth errors by less than 35% and decreases 
from 221 to 137. We expect to be able to determine 
the calibration to even better than 1% from comparison 
of ground-based data with high-resolution space-based 
images over a small fraction of the total survey area. 

The imaging data from which shear maps are derived 
can also be used to infer galaxy power spectra. The large- 
scale feature from matter-radiation equality used here 
and the baryonic oscillations at smaller scales can also 
be used to infer distances miii. 

The distance-redshift and growth-redshift relations 


provide two observational windows on the physics of ac¬ 
celeration. While we have illustrated the utility of a sec¬ 
ond window with a specihc example, the extra informa¬ 
tion may prove crucial to the unraveling of the mystery 
of acceleration in ways we have not yet imagined. 
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